Brillouin Scattering Study of Propylene Carbonate: An Evaluation of 
Phenomenological and Mode Coupling Analyses 
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Brillouin scattering spectra of the molecular glassformer propylene carbonate (PC) in the tem- 
perature range 140 K to 350 K were analyzed using both the phenomenological Cole-Davidson 
memory function and a hybrid memory function consisting of the Cole-Davidson function plus a 
power-law term representing the critical decay part of the fast /3 relaxation. The spectra were also 
analyzed using the extended two-correlator schematic MCT model recently employed by Gotze and 
Voigtmann to analyze depolarized light backscattering, dielectric, and neutron-scattering spectra 
of PC [W. Gotze and Th. Voigtman, Phys. Rev. E 61, 4133 (2000)]. We assess the ability of the 
phenomenological and MCT fits, each with three free fitting parameters, to simultaneously describe 
the spectra and give reasonable values for the a-relaxation time r a . 

PACS numbers: 64.70.Pf, 78.35.-fc, 67.40.Fd 



I. INTRODUCTION 

Brillouin scattering spectroscopy has been one of the 
experimental techniques widely applied to the study of 
glass-forming materials in the continuing effort to achieve 
a complete understanding of the liquid-glass transition. 
The evolution of the Brillouin spectrum with tempera- 
ture reflects the interaction of longitudinal sound waves 
(q ~ 10 5 cm -1 ) with structural relaxation. Typically, as 
T decreases from above the melting temperature Tm to 
below the glass-transition temperature Tg, the Brillouin 
linewidth Aw^ first increases, passes through a maxi- 
mum, and then decreases again, while the Brillouin peak 
position lob shifts monotonically to higher frequencies, 
with the increase occurring predominantly in the tem- 
perature range where Au>b is biggest. 

In the earliest experiments, only the temperature de- 
pendence of the Brillouin shift and linewidth were stud- 
ied. Subsequently, the full spectrum was analyzed and 
compared with theoretical predictions derived primar- 
ily from generalized hydrodynamics models. Either the 
Mountain formalism [Mountain 1968] or a generalized 
damped harmonic oscillator model were used in which 
the memory function, i.e. the frequency-dependent part 
of the longitudinal viscosity, represents the structural re- 
laxation dynamics. Generally, acceptable fits could be 
obtained using an exponential or stretched-exponential 
empirical memory function m(t) (or a Cole-Davidson 
memory function ro[w]), with the relaxation time t{T) 
treated as an adjustable fitting parameter. However, 
the t(T) values determined from such fits were usually 
found to increase much more slowly with decreasing T at 
low temperatures than the values determined with other 
techniques such as dielectric spectroscopy or ultrasonics. 

The origin of this apparent disagreement can be at- 
tributed to the use of these "a-relaxation only" models 



for m(t) which ignore the fast part of the relaxation pro- 
cess preceding the final a-relaxation. The need for an 
"additional fast contribution" to m(u>) was first noted by 
Loheider et al [Loheider 1990] (also see [Cummins 1994]). 
Several empirical approaches to extending m{oS) or m(t) 
to include the missing fast j3 relaxation have been de- 
scribed in the literature, including: (1) adding a sec- 
ond Debye, Cole-Cole, or Cole-Davidson function to 
the Cole-Davidson a-relaxation term [Soltwisch 1998] 
[Monaco 2000, Monaco 2001]; (2) extracting a memory 
function including both a and (3 relaxation from depo- 
larized backscattering spectra and using it for fitting the 
Brillouin spectra [Cummins 1994] (a procedure that suc- 
ceeded for CKN but not for other materials); (3) tak- 
ing m(i) as the sum of a stretched exponential plus a 
damped high-frequency oscillatory term [Lebon 1997]; 
and (4) adding a term Buj a to the Cole-Davidson func- 
tion to approximate the critical-decay component of the 
fast /3-relaxation predicted by the mode coupling theory 
(MCT) [Du 1994, 1996, Cummins 1994, Pick 2000]. This 
additive "hybrid model" approach, initially suggested by 
Gotze [Gotze 1993], has recently been shown to pro- 
vide excellent fits to depolarized backscattering spectra 
of toluene [Wiedersich 2000] and has also been used in the 
analysis of depolarized Brillouin spectra of metatoluidinc 
[Pick 2000]. 

With the additional fast relaxation included through 
any of these approaches, the a-relaxation time T a is found 
to increase much more rapidly with decreasing temper- 
ature than it does with a-only models, consistent with 
other experimental measurements. However, while this 
approach can resolve the t{T) disagreement, it usually 
requires additional fitting parameters and does not ex- 
plain the origin of the structural relaxation dynamics 
that it describes. Furthermore, various authors have cho- 
sen quite different empirical fitting functions and have 



usually obtained good fits to Brillouin spectra, suggesting 
that this procedure is not able to distinguish effectively 
between different models. 

Another approach to the analysis of Brillouin scatter- 
ing spectra utilizes the mode coupling theory (MCT) in 
which the slowing down of structural relaxation with de- 
creasing T arises naturally from nonlinear interactions 
among density fluctuations. A major advantage of this 
approach is that both the a and (3 regions of the relax- 
ation dynamics occur spontaneously in the solutions of 
the MCT equations. However, while MCT can provide 
a quantitative microscopic description of liquid dynam- 
ics in principle, it requires the intermolccular potential 
[or the static structure factors S(q)\ as input, and this 
information is not generally available for most real glass- 
forming materials. Therefore, MCT analyses are usually 
restricted to comparison with asymptotic predictions in 
which the detailed microscopic interaction information 
is absorbed into a single parameter, the exponent pa- 
rameter A. Alternatively, a simplified "schematic MCT" 
can be used in which the full set of density correlators 
is replaced by a small number (usually one or two), and 
the coupling constants are treated as adjustable fitting 
parameters. This schematic approach captures the ba- 
sic mechanism underlying the MCT more fully than the 
asymptotic analyses, but sacrifices the material-specific 
microscopic interpretation of the parameters. 

A convenient schematic MCT model, initially intro- 
duced by Sjogren [Sjogren 1986], uses two correlators: 
</>(£) for the "system" and <fi s (t) for the "probe", i.e. 
the variable being measured in a particular experiment. 
In this Sjogren model, the interesting critical behav- 
ior occurs in <j>(t) while the temperature dependence 
of 4> s (t) results from its being coupled to 4>(t). This 
model was successfully used to analyze depolarized light- 
scattering spectra of glycerol [Franosch 1997] and OTP 
[Singh 1998]. Ruffle et al [Ruffle 1999] showed that Bril- 
louin scattering spectra of NaLi(P03)2 could also be an- 
alyzed with the Sjogren model, using parameters for the 
system correlator <p(t) determined from inelastic neutron 
scattering spectra. 

Recently, Gotze and Voigtmann have extended the 
Sjogren model by including hopping parameters A and 
A s in the equations of motion for both <fi(t) and 4> s (t), 
which allows the analysis to be extended to low tempera- 
tures where relaxation is dominated by activated hopping 
processes. They analyzed dielectric, depolarized light 
scattering, and neutron scattering data for propylene car- 
bonate (PC), using a single set of system parameters for 
4>(t), optimized simultaneously for all the experiments 
[Gotze 2000] . This result suggests the possibility of an- 
alyzing PC Brillouin scattering spectra using the system 
parameters already determined by Gotze and Voigtmann. 
Moreover, if the probe hopping parameter A s , frequency 
tt s , and damping v s are taken as fixed by the depolar- 
ized light scattering fits, then the analysis will include 



only one adjustable MCT fitting parameter, the coupling 
constant V s , which determines the strength of coupling 
between 4>{t) and (j> s (t). The fits should then simultane- 
ously describe the Brillouin spectra and also determine 
the value of the structural relaxation time r Q . Such an 
analysis could provide a stringent test of the utility of 
MCT for the analysis of Brillouin scattering spectra. A 
test of this procedure was the principal motivation for 
the present Brillouin scattering study of PC. 

We describe the experiments in Sec. II. We analyze the 
data in Sees. Ill and IV using two empirical memory- 
function models, first a conventional a-relaxation-only 
Cole-Davidson (CD.) function and then the hybrid 
model. (A preliminary version of the CD. and hybrid 
fits was presented by Frank [Frank 2000].) In Sec. IV we 
also include a discussion of the role of thermal diffusion, 
followed by a second hybrid model analysis in which the 
CD. stretching parameter j3 is constrained to be equal 
to the MCT critical exponent b. In Sec. V we present our 
analysis of the spectra with the two-correlator extended 
schematic MCT approach. In Sec. VI, we provide discus- 
sion and conclusions including a brief evaluation of the 
ability of these and related approaches to extract mean- 
ingful information about structural relaxation dynamics 
from the analysis of Brillouin scattering spectra. 



II. EXPERIMENT 



A. Material 



Propylene carbonate (PC; Tq = 160 K) is a molecular 
glassforming material that has been extensively studied 
with a wide range of experimental techniques including 
Brillouin scattering [Elmroth 1992; Du 1994]. In Table I 
we list the principal properties of PC that are relevant 
for the present study. 

A recent dielectric spectroscopy study of PC by Schnei- 
der et al [Schneider 1999] determined the complex dielec- 
tric constant e(u>) over a frequency range of 18 decades 
for temperatures from 153 K to 293 K. Cole-Davidson fits 
to the a peaks in the e"(u>) spectra produced the t^(T) 
values shown in Fig. 1 by the open circles. Schneider et 
al carried out fits to these values with several empirical 
fitting functions. Their extended free volume fit is shown 
in Fig. 1 by the solid line. 

Du et al [Du 1994] analyzed the a peaks in x"( w ) spec- 
tra, obtained from their depolarized PC backscattering 
spectra, with a Kohlrausch function. The resulting t^ l 
values, shown in Fig. 1 by squares, have temperature de- 
pendence similar to the dielectric t^{T) although they 
are approximately three times smaller. 
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FIG. 1. (TAU). Propylene carbonate (PC) relaxation 
time r a (T). Open circles and solid line: dielectric 
measurements and extended free-volume fits [Schnei- 
der 1999]. Squares: depolarized light scattering spectra 
[Du 1994]. Solid circles: Brillouin scattering t^ d (T) 
from Cole-Davidson fits (Sec. III). Broken line: di- 
electric free volume fit multiplied by 0.085 to match 
Brillouin T a (T) in the region of the maximum Bril- 
louin linewidth. *: t^ ict obtained from the extended 
schematic MCT fits. 



The analysis of our Brillouin scattering spectra using 
a Cole-Davidson memory function, described in the next 
section, resulted in the t£ d values shown by the solid cir- 
cles in Fig. 1. As mentioned in the Introduction, at low 
temperature these values do not increase as rapidly with 
decreasing temperature as the dielectric and backscatter- 
ing results, suggesting that the Cole-Davidson function 
provides an incomplete representation of structural re- 
laxation. (At T — Tq, t^ d has only reached ~ 30 ns, 
while typically T a (Tc) is ~ 100 s.) 

In our 90° PC Brillouin scattering experiment, the 
maximum Brillouin linewidth occurs at T ~ 230 K. At 
that temperature the peak of the a relaxation coincides 
with the Brillouin line, so wbt„ ~ 1 which gives t^ a 
(230) ~ 3.2 x 10~ n sec, a value ~ 12 times smaller 
than the t^ of Schneider et al. The difference presum- 
ably arises from the fact that t£ a characterizes longitu- 



dinal displacement dynamics while t^ characterizes ori- 
entational dynamics. In Fig. 1, the lower line is the free 
volume fit of Schneider et al divided by a factor of 11.75 
to match t„ in this region: 



logio WT) = -A + B/{T - To + [(T - T Q y 



CT] 1 ' 2 } 
(1) 



with A = 12.56, B = 309, C = 4.82, and T = 162 K 
(t„ in seconds). In constructing memory functions for 
data analysis beyond the Cole-Davidson model, we will 
use Eq. (1) to fix the a-relaxation time of the longitudinal 
acoustic mode t^ a . 



B. Sample Preparation 



Propylene carbonate (anhydrous, 99.7 %) was pur- 
chased from Sigma- Aldrich. PC was loaded in a glove box 
under dry nitrogen atmosphere into a distillation flask fit- 
ted with a stopcock. The flask was then transferred to a 
vacuum distillation system and distilled at ~ 110 C into 
glass sample cells that were flame sealed under vacuum. 
A sample cell was installed in an Oxford LN2 coldfinger 
cryostat with an ITC-4 temperature controller for the 
Brillouin scattering experiments. 



C. Brillouin scattering experiments 



Spectra were collected in both 90° VV geometry and 
174° VH near-backscattering geometry with a Sander- 
cock 6-pass tandem Fabry-Perot interferometer at tem- 
peratures from 140 K to 350 K in 5 K or 10 K steps. 
The 514.5 nm single-mode argon laser power at the sam- 
ple was typically 170 mW. The interferometer finesse was 
~ 90; its contrast was better than 10 7 . Complete sets of 
spectra were collected with mirror separations d = 2, 5, 
and 10 mm for VV and d — 2 mm and 10 mm for VH. 
Spectra were accumulated during 1 to 2 hour runs. For 
the d = 10 mm separation spectra two 2-hour runs were 
collected and averaged to improve the signal-to-noise ra- 
tio. In Fig. 2 we plot the Stokes side of the three VV 
spectra (a) and the two VH spectra (b) for T = 220 K. 
(The dark counts have already been subtracted.) 

In order to test for possible experimental artifacts, 
Iiso(u) PC Brillouin spectra at several temperatures 
were collected at both the Technical University of Mu- 
nich and at the City College of New York with different 
samples, different tandem Fabry-Perot spectrometers, 
and different experimental procedures. The spectra 
obtained were in very good agreement. The spectra 
shown in this work were all collected in New York. 
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FIG. 2. (EXP1). (a) Stokes sides of 220 K PC 90° VV 
Brillouin spectra with d = 2, 5, and 10 mm. (b) 174° 
VH spectra with d — 2 mm and 10 mm. (c) The com- 
posite 90° VV spectrum obtained by shifting the 2 mm 
and 5 mm spectra of (a) to optimize overlap with the 
10 mm spectrum, the composite VH spectrum obtained 
similarly and then shifted to optimize overlap with the 
VV spectrum for frequencies above 30 GHz, and the 
difference spectrum Iiso(j-^) obtained by subtraction. 
The background from dark counts has already been 
subtracted in the spectra of (a) and (b). 



D. Data reduction: extraction of the density 
fluctuation spectra Iiso(u) 



The VV polarized Brillouin spectra contain contribu- 
tions from density fluctuations [Iiso(oj)] an d also from 
oricntational dynamics and collision-induced scattering 
[Iani(u)]- Depolarized (VH) backscattering spectra are 
due to Iani{w) only, and can therefore be used to remove 
the Iani{w) contribution from the VV spectra. 

Fioretto et al [Fioretto 1999] [Monaco 2001] have car- 
ried out the subtraction by combining spectra obtained 
with different mirror separations to obtain broadband 
VV and VH spectra, and scaling the VH spectra to over- 
lap the VV spectra for frequencies above ~ 40 GHz where 
IiSo(w) should be negligible. Subtracting the scaled 



Ivh(w) from the Iyy(u)) then provides Iisoi 1 ^)- We have 
also followed this procedure. 

The Stokes sides of the d = 2 and 10 mm 220 K VV 
spectra, shown in Fig. 2(a), were rescaled vertically to 
optimize overlap with the d — 5 mm GHz spectrum, 
producing the composite 90° Iyy{ui) spectrum shown in 
Fig. 2(c). The 9 = 184° Ivh{u) spectra were similarly 
combined, and the composite Ivh spectrum was again 
scaled vertically to match the Iyy spectrum in the 30 - 
50 GHz region. Finally, the scaled Iyjj spectrum was 
subtracted from the Iyy spectrum to produce the differ- 
ence spectrum Iiso(^>) as shown in Fig. 2(c). 

The full set of Iiso(w) spectra for temperatures from 
140 K to 350 K obtained with this procedure is shown 
in Fig. 3. For each set of VV spectra, we also recorded 
an instrument profile spectrum with d = 10 mm. The 
width of the instrumental profile was typically FWHH 
= 0.17 GHz. To reduce additional Brillouin peak broad- 
ening due to the finite input aperture (which leads to 
collecting light scattered at slightly different angles), the 
input aperture was reduced sufficiently to decrease the 
additional broadening at the Brillouin frequency to less 
than the instrumental width. This additional geometrical 
broadening (typically < 0.13 GHz) was then calculated 
for every temperature from the known collection geome- 
try and convoluted into the instrument function. In the 
fitting procedures described in this paper, the instrument 
function so obtained was convoluted with the theoreti- 
cal spectrum for comparison with the Iiso(u) spectrum. 
The fits were carried out with a conventional nonlinear 
least-squares fitting program (NLLSQ). 

We first fit each Brillouin peak in the d = 10 mm 
VV spectra to the damped harmonic oscillator func- 
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FIG. 3. (EXP2). Complete set of composite Iiso{w) 

spectra for 28 temperatures from 140 K to 350 K ob- 
tained by the procedure shown in Fig. 2. T = 140, 150, 
160, 165, 170, 180, 185, 195, 200, 205, 210, 215, 220, 
225, 230, 235, 240, 250, 260, 270, 280, 290, 300, 310, 
320, 330, 340, and 350 K, from right to left. 



J( W ) = 



— Xm[tu — tu B — ioj^b] 
to 

_ hlB 



(2) 



including convolution with the instrument function. In 
Fig. 4 the resulting values of wb(o) and 7b (6) are shown 
by the points. For a simple liquid, both lob and 73 
would have simple monotonic temperature dependence. 
The observed temperature dependence of o>b and 7b , as 
noted in the introduction, is the signature of structural 
relaxation dynamics as seen in Brillouin scattering. The 
maximum in 75 at T ~ 230 K indicates that the a peak 
moves through the Brillouin line at that temperature. 



III. DATA ANALYSIS 1: THE COLE-DAVIDSON 
MODEL 

A. The memory function approach 

The isotropic Brillouin spectrum Iiso (w) is usually at- 
tributed to density fluctuations p q (t) where the scattering 
vector q = (2wn/X) sin(0/2). The spectrum consists of a 
triplet: the central thermal diffusion mode (which is gen- 
erally located at frequencies too low to observe in these 
experiments and will be ignored in most of the fits) and 
the Brillouin doublet. The spectrum is conventionally 
described by 



I(w) 



Jb[7o+ro"(w)] 



[lu 2 — Wq + wm'(w)] 2 + [lujq + wto"((jj)] 5 



(3) 



[The particular Laplace transform convention used to de- 
rive Eq. (3) is discussed in the Appendix.] 

Equation (3) is very general; m(w) = m'(u>) + im"(u)) 
is the Laplace transform of the memory function m(t) 
which can be constructed to include structural relax- 
ation, thermal diffusion, and translation-orientation cou- 
pling for liquids of anisotropic molecules [Chung 1971; 
Boon 1980; Scopigno 2000]. Usually, however, m(t) 
represents the relaxing longitudinal viscosity, hrst in- 
troduced by Mountain [1968] to describe the damping 
of sound waves by interaction with internal degrees of 
freedom of the molecules, and later utilized to represent 
the interaction of sound waves with structural relaxation. 
(Both mechanisms may contribute, however, as discussed 
below.) In the long-wavelength limit, cjo = C$q where 
Co is the limiting low-frequency adiabatic sound veloc- 
ity, and 70 is the "regular" sound attenuation coefficient. 
Note that Eq. (3) is the equation for a damped harmonic 
oscillator with a frequency-dependent damping function 
j(u>) = 70 +m(w). 

Ideally, the Iiso(<*>) spectra should be analyzed to 
yield m{uj) directly. But in practice this is not possible 
because the spectra only cover a range of two decades 
at most, while m(u>) extends over many decades. Also, 



the spectra are dominated by the Brillouin components 
[which would be the complete spectrum of Eq. (3) if 
m(u>) = 0], and the modification of I(u>) produced by in- 
clusion of m(w) is not very sensitive to the detailed form 
of m(u>). Therefore, data analysis almost always proceeds 
by selecting a parameterized model for m{u>) and varying 
the parameters at each temperature to optimize the fits. 

The primary requirements on m(u>) are that m'(tUB), 
the real part of m(uj) at the frequency of the Brillouin 
peak, must be adjusted to shift the Brillouin peak from 
luq to u>b; and m"(u>B), the imaginary part of m(u>) at 
the frequency of the Brillouin peak, must be adjusted so 
that 7o+m"(tc>B) will give the correct Brillouin linewidth. 
If these two conditions are simultaneously satisfied, then 
the particular form chosen for m(oj) will only influence 
the details of the lineshape. Therefore the essential re- 
quirement is that m" (u>b) I m' (ujb) must have the "cor- 
rect" value. Nevertheless, this is not a trivial requirement 
since m'(uj) and m"(uj) are connected by Kramcrs-Kronig 
relations. 

The prototype memory function is the exponential 
e~*' T , multiplied by a coupling constant A 2 . (It is usu- 
ally called Debye relaxation although it was first used by 
Maxwell in his theory of viscoelasticity.) The simplest 
memory function m(t) for structural relaxation that fits 
various experimental data in the frequency region of a 
relaxation is the stretched-exponential (KWW) function 
e -(*/ r ) K _ An approximation to the Fourier transform of 
this m(t) is provided by the Cole-Davidson (CD.) func- 
tion 



l(«) = {A 2 /u;)[{1-iu;t)- (j - 1] 



(4) 



which has frequently been used, together with Eq. (3), to 
analyze Brillouin-scattering spectra. In practice, the fit- 
ting procedure adjusts the coupling constant A 2 so that 
the real part of m(w) produces the correct lob for the 
Brillouin peak, while r is adjusted to give the correct 
Brillouin linewidth. 

There are several problems that should be noted in 
carrying out this analysis, as well as analyses with more 
elaborate memory functions. 

(1) Many glassforming materials are molecular liquids, 
and internal vibrational modes may also contribute to 
m(uj), as originally proposed by Mountain [1968]. Re- 
cently, Monaco et al ([Monaco 2000]) have carried out a 
Brillouin scattering study of OTP from which they con- 
clude, for OTP, that the fast part of the relaxation pro- 
cess is entirely due to intramolecular vibrational modes. 
The assignment results from the presence of additional 
low frequency structure (Mountain mode) in VV, but 
not VH, spectra, in both the glass and the crystal. This 
conclusion has recently found additional support from a 
molecular dynamics study of OTP [Mossa 2001]. There 
is no evidence, however, that it applies to other materials 
such as PC, and it will not be considered here. 



(2) The identification of Iiso{w) with density fluc- 
tuation rests on the assumption that scattering due to 
orientational fluctuations, which is typical for molecular 
glassforming materials, is eliminated by the subtraction 
of I an i(oj) from Iyy{uS) described in the previous section. 
However, two consequences of anisotropy may remain af- 
ter subtraction. First, since longitudinal current involves 
both compression and shear, rotation-translation cou- 
pling modifies the longitudinal current correlator [which 
determines the spectrum of p q (t)] as shown in Eq. (A. 14) 
of Ref. [Dreyfus 1999]. Second, as noted recently by Latz 
and Letz [Latz 2000], because of this coupling the ori- 
entational dynamics also reflect the longitudinal acoustic 
mode, and this part of the orientational dynamics is not 
removed by subtraction. Therefore, both density fluctu- 
ations and orientational dynamics may contribute to the 
Iiso{w) spectra. 

While the extent of these corrections due to molecular 
anisotropy is not yet clear, preliminary simulations, car- 
ried out in collaboration with R.M. Pick, indicate that 
they are relatively minor. We will therefore use the 
conventional memory function formalism in this paper, 
recognizing that the values of the parameters used to 
describe the frequency-dependent longitudinal viscosity 
may reflect some aspects of translation-rotation coupling 
and molecular anisotropy and may also include some con- 
tribution from intramolecular vibrational dynamics. 



B. Cole-Davidson Fits 



We first analyzed the Iiso (w) spectra of Fig. 3 using 
Eq. (3) with the CD. memory function m(ui) of Eq. (4). 
The fit results, for a subset of the spectra, are shown in 
Fig. 5 [note (IIIB.l)]. As often found in the past, excel- 
lent fits can be obtained if all the parameters are kept 
free. However, it is preferable to fix as many parame- 
ters as possible using independently determined values, 
since the parameters are rather strongly coupled. (This 
coupling can also lead to unstable parameter values in 
the fit results.) Here we consider each of the relevant 
parameters separately. 

(1) (3. The fits are relatively insensitive to this param- 
eter. We took (3 — 0.68 (corresponding to (3r — 0.77) as 
found in the analysis of depolarized backscattering spec- 
tra by Du et al [Du 1994]. 

(2) u>o(T). Since luq = C$q where Cq is the limit- 
ing low-frequency adiabatic sound velocity and q is the 
scattering wavevector (Awn/X) sin(0/2), u)q(T) can be 
fixed if Cq(T) and the refractive index n(T) are known. 
Co(T) was previously measured by ultrasonic experi- 
ments at 5 MHz and 15 MHz by Du et al [Du 1994] 
which gave C (T) = 2.507 x 10 5 - 361.2 T cm/sec. 



| Wo. 1 May | 



1 1 - 


• • .. 








(a) 


10- 


^ quadratic 


•. 








9 - 












X 

a 

B 8- 
CM 


linear" 1 "** 


\. *•. 








S 

7 - 






^- 






6 - 












5 - 


TG 

1 


TM 








I 
150 


I 
200 


I 
250 
T(K) 


300 


I 

350 


1 .2- 




_• 


•• 




(b) 


1 .0- 




' 


• 






■ST 0.8- 

I 

O 

& 0.6- 
0.4- 








• 


" • • . 


0.2-. 


. ..'" 
























150 


200 


250 


300 


350 






T(K) 







FIG. 4. (Wo). Fits of the d = 10 mm 90° VV Brillouin 
peaks to Eq. (2) including convolution with the instru- 
ment function, (a) wb/27t(GHz) (circles). Short solid 
line: coq from ultrasonic and refractive index data; long 
solid line: loq increased by 3.5% to match lob in the 
high-temperature region. Solid curved line: quadratic 
fitw (r) = 13.211 -3.291 x 10~ 2 T +2.6973 xlO~ 5 T 2 . 
(b) 7b /27t (GHz) (circles). Solid line: 7o = T/1000. 

However, the ultrasonic data only cover the range 240 
to 293 K. At lower temperatures, there is significant dis- 
persion in the ultrasonic sound speed, while higher tem- 
perature measurements are technically difficult. The re- 
fractive index n(T) = 1.5314 - 3.752 x 10" 4 T (see Ta- 
ble I) corresponds to the sodium D line and may be dif- 
ferent at 514.5 nm. The ojq(T) values implied by these 
measurements are shown in the upper panel of Fig. 4, by 
the short solid line. A linear extrapolation of these val- 
ues, shifted up by 3.5 %, is indicated by the long solid line 
labeled "linear": uj (T)/2t: = (10.71 - 0.0165 x T) GHz. 
The shift, which is within the experimental errors of the 
Cq(T) and n(T) data, was chosen so that at high temper- 
atures the Brillouin peak frequency lub is equal to luq, as 
expected. Note, however, that the high-T lob values sug- 
gest some small upward curvature away from the linear 
extrapolation above 320 K. 

Initially, fits were carried out with this linear- T esti- 
mate of u)q(T). However, the fits were not all satisfactory, 



in part because of the high-T curvature in lob seen in 
Fig. 4. For materials where Cq(T) measurements extend 
over larger temperature ranges, there can be significant 
upward curvature to Cq(T). We therefore tried fitting 
u>o(T) to the six highest temperature u>b(T) values plus 
one u>o(T) value obtained by a free CD. fit at 230 K 
where the Brillouin linewidth is largest. The result was 
the quadratic u>o(T) function shown by the upper solid 
line in Fig. 4 labeled "quadratic" : 



w (T) = 13.211 - 3.291 x 1CT 2 T 



2.6973 x 1CT 5 T 2 
(5) 



which was used in all subsequent fits. 

(3) 7o(r). This term is conventionally included to 
represent "regular" damping of the sound waves by an- 
harmonic processes not related to structural relaxation. 
It can also represent a first frequency-independent ap- 
proximation for the fast part of the structural relax- 
ation not included in a-relaxation-only models such as 
the CD. function. Usually, 70 is taken as a temperature- 
independent constant, fixed from the Brillouin linewidth 
at low temperatures. 

However, since anharmonic damping processes often 
tend to become stronger with increasing temperature, 70 
may increase with T. Also, at temperatures well above 
the temperature of the maximum Brillouin linewidth, the 
linewidth due to 70 plus structural relaxation is approx- 
imately given by Aujb ~ 7o + A 2 r/3. Fits with constant 
70 often give A 2 values that increase at high tempera- 
tures, a result that appears unphysical. This apparent 
increase in A 2 can be avoided by allowing 70 to increase 
with increasing T. To avoid introducing additional fit- 
ting parameters we arbitrarily assumed 70 to be a linear 
function of T, and took 70 (T) = T/1000 as shown by the 
solid line in the lower panel of Fig. 4. 

(4) A 2 (T). Since at high temperatures u>b — Wo while 
Aujb = 7o + A 2 r/3, the effective fitting parameter for 
Aujb is the product A 2 r. So if A 2 and r are both free 
parameters, the fits tend to be unstable. We therefore 
carried out the fits for temperatures up to 235 K with 
the three free fitting parameters: A 2 , r, and the scale 
factor Iq. For T > 235 K, A 2 was fixed at its value 
at ~ 235 K, and the fits were carried out with only r 
and Iq free. The parameters t(T) and A 2 (T) obtained 
from these fits are shown in Table II along with the re- 
duced x 2 value for each fit. The CD. fits, shown in 
Fig. 5, are satisfactory although better fits could be ob- 
tained by keeping more parameters free, especially ojq. 




FIG. 5. (CDFit). Cole-Davidson fits to the Iiso{u) 
spectra of Fig. 3 to Eqs. (3) and (4), with (3 = 0.68, 
lu /2it = 13.211 - 3.291 x 10~ 2 T +2.6973 x 10~ 5 T 2 , 
7o/27r = T/1000. The fits are shown for T = 150, 160, 
170, 180, 195, 210, 220, 230, 240, 260, 280, 300, 320, 
and 340 K. 

However, the major point of the CD. fits, as shown 
by the solid circles in Fig. 1, is that for temperatures be- 
low ~ 230 K, t£ d determined from these fits does not 
increase very rapidly with decreasing T. At T = Tq = 
160 K, t^' d found from the CD. fits is only ~ 30 ns, com- 
pared tor~ 100 s found in other experiments. This dis- 
agreement is generally recognized as indicating that the 
CD. model (or any a-relaxation-only model) for m(u)) 
is incomplete. Once the a peak moves well below the 
frequency of the Brillouin peak, m"(uj) of the CD. func- 
tion [with reasonable t(T) values] decreases rapidly with 
decreasing temperature, while the shape of the spectrum 
indicates that a significant contribution from m"{uS) is 
still present in the region of the Brillouin peak. 

The fact that the CD. function is an incomplete repre- 
sentation of the structural relaxation dynamics can also 
be seen directly in the PC dielectric data of Schneider 
et al [Schneider 1999]. Their Fig. 2 shows CD. fits to 
the dielectric data which are excellent for the a peak, 
but which fall below the experimental data in the high- 
frequency wing. 



IV. DATA ANALYSIS 2: INCLUDING THE FAST 
(3 RELAXATION 

The t{T) values obtained from the Cole-Davidson fits 
exhibit temperature dependence at low temperatures 
that disagrees strongly with the results of other exper- 
iments, as seen in Fig. 1. As mentioned in the intro- 
duction, this disagreement is a characteristic failure of 
a-relaxation-only models for m(u)), which indicates the 



need for extending m(u>) to include fast relaxation pro- 
cesses. Several authors have proposed models that ex- 
tend the CD. model empirically by adding another phe- 
nomenological term to Eq. (4), e.g. a Cole-Cole term 
[Soltwisch 1998], a damped oscillator [Lcbon 1997], or a 
Lorentzian [Monaco 2000B]; however, attributing a fast 
relaxation time 173 to structural relaxation is incompat- 
ible with MCT which describes the fast relaxation by 
power laws in u) or t for which there is no 773 [Gotze 1992] . 
A major advantage of this approach is that the additional 
fast CD. or Lorentzian, with relaxation time 773 ~ 30 ps, 
is able to explain the extra "Mountain mode" seen in 
OTP at low temperatures [Monaco 1999, Fig. 1(a)]. 
(Recently, however, this feature has been attributed to 
Mountain's original mechanism of intramolecular dynam- 
ics rather than to structural relaxation [Monaco 2000, 
Mossa 2001].) 



A. The hybrid model 

We have utilized a different additive memory function, 
the hybrid model (mentioned in the introduction), which 
has been employed in several previous Brillouin scatter- 
ing studies. We begin with the CD. function [Eq. (4)] 
and add a term to m(uj) proportional to uj a ~ l which rep- 
resents the t~ a "critical decay" part of the fast (3 process 
in MCT. This superposition approximates the two-step 
relaxation scenario of MCT, and does not introduce an- 
other relaxation time 773. 

In some previous studies, we have used a related proce- 
dure to analyze depolarized backscattering spectra. For 
orthotcrphenyl, we combined extended MCT hts for the 
/3-relaxation region with KWW fits for the a peak [Cum- 
mins 98] . A comparison of these two hybrid procedures 
has not yet been attempted. 

Two problems arise when the term t~ a is added to 
m(t). First, the resulting u> a term in ujm{uj) extends to 
arbitrarily high frequencies. This is presumably unim- 
portant for Brillouin scattering, however, since fits to the 
Brillouin spectra are limited to u>/2ir < 15 GHz. Sec- 
ond, m"(w), which appears in the numerator of I(u>) in 
Eq. (2), diverges at low frequencies as u a ~ l resulting in 
the introduction of a spurious "central peak" . To elimi- 
nate this artifact, we have included an exponential cutoff 
in the critical decay part of m(t) [Gotze 2000B]: 

-(t/r) 



m crit (t) = Bt a e 



(6) 



where r is the same t„ of the CD. function, Eq. (4). 
(Note that this cutoff factor has not been included in 
previous analyses based on the hybrid model.) 

With this modification, the hybrid memory function 
becomes 



[see the Appendix for the derivation of Eq. (7) .] 

At high frequencies where lut >> 1, the critical decay 
part of Eq. (7) becomes 



ujm C R{uj) = [Bi + iB 2 ]oJ a 



(8) 



{where B x = Bsin[7r(a - 1)/2]T(1 - a) and B 2 = 
Bcos[ir(a — 1)/2]T(1 — a) }. Equation (8) is the familiar 
high-frequency power-law behavior predicted by MCT. 

In Fig. 6 we plot utm"(uj) of Eq. (7) with r = 0.3 ns, 
(3 = 0.68, and a = 0.29. The broken line is the 
Cole-Davidson term only, and the lower solid curve is the 
critical decay term only, including the exponential cutoff. 
The upper solid curve is the full hybrid memory function 
wm"(w) of Eq. (7) with B/A 2 = 0.005. This figure illus- 
trates how the resulting wm"(w) goes over smoothly from 
the CD. a-peak behavior at low frequencies to the ut a be- 
havior at high frequencies, with a minimum between the 
two regions. The inset to Fig. 6 shows the full m"(oj) 
with and without the e~ t ' T cutoff, illustrating how the 
spurious central peak is eliminated by the exponential 
cutoff. 



B. The critical exponent a 



In MCT, the function Lom"(uj) is predicted to have 
minimum between the high-frequency von-Schweidler 
wing of the a peak and the critical decay region simi- 
lar to that in Fig. 6. In the region of this minimum, the 
MCT interpolation equation gives u>m"(u>) in terms of 
the critical exponents a and b by 



um"(u) = Lum l } nin (uj)[b{uj/w m in) a 
+a(u mm /uj) b }/(a + b) 



(9) 



wm(w) = A 2 [(l - iojT)- 13 - 1] 
+iujBT{l-a){T- 1 -iu) a - 1 



(7) 



where a is the critical exponent and b is the von- 
Schwcidler exponent, representing the high-frequency 
wing of the a peak. 

Since in our hybrid memory function, Eq. (7), the CD. 
function for lot >> 1 is oc oj~P, the interpolation Eq. (9) 
with b = [3 will determine an "effective" a e ff to use in 
Eq. (7) at temperatures high enough for the minimum to 
be observable (from previous MCT fits, b w 0.50 while 
[3 ~ 0.68). At lower T, where the minimum has disap- 
peared from the frequency region of the Brillouin spec- 
trum, uim"{uj) may have the form of a power law in to, but 
since this region no longer corresponds to the asymptotic 
region of the minimum, the apparent power law behavior 
of uim"{u)) oc Lu a "ff can produce an a e // quite different 
from the critical exponent a itself. 
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FIG. 6. (HYB1). The hybrid memory function wm"(w) 
of Eq. (7) with r = 0.3 ns, = 0.68, and a = 0.29. The 
broken line is the CD. function and the lower solid line 
is the critical decay term. The ratio of the critical decay 
to CD. contributions is B/A 2 = 0.005. Inset: m"(w) 
vs lu for the hybrid model [Eq. (7)] without a cutoff 
(solid line), with r = 0.01 ns, A = 20 GHz, and B = 5; 
m"{uj) with the cutoff e~'/ T in Eq. (6) (broken line). 
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FIG. 7. (HYBFIT2). VH backscattering spectra mul- 
tiplied by frequency with fits to the MCT interpolation 
equation [Eq. (9)] (using b = /3 = 0.68) (T = 170, 180, 
185 K) or the power law uj a 'ff (T = 140, 150, 160 K). 
The values of a e ff determined in these fits were used in 
the low-temperature Iiso{u) hybrid model fits shown 
in Figs 8 and 9. 



We utilized the depolarized backscattering spectra 
Ivh(w), multiplied by lu to estimate u>m"(u)) in order 
to find a e ff(T). In Fig. 7 we show six u>Ivh{w) spectra 
for T = 140, 150, 160, 170, 180, and 185 K. The spectra 
for T — 170, 180, and 185 K, where the minimum is vis- 
ible, have been fit to the interpolation equation [Eq. (9)] 
using /3 = 0.68 for b, while the spectra for T = 140, 150, 
and 160 K have been fit to Lo a "ff . The values of a e ff(T) 
found from these fits are: 



T 


a eff 


.40 


0.62 


.50 


0.59 


.60 


0.52 


.70 


0.42 


.80 


0.30 



In the fits to our Brillouin spectra, we used 
these values of a e ff(T) for temperatures up to 
180 K and a = 0.29 for all higher temperatures. 



C. Minimum free parameter analysis 

We initially tried to fit the spectra with Eq. (3) and the 
hybrid memory function Eq. (7) with the smallest pos- 
sible number of free parameters. We fixed loq{T) with 
Eq. (5), fixed 70 (T) = T/1000 as in the Cole-Davidson 
fits, fixed r using Eq. (1), fixed /3 = 0.68, and fixed the 
ratio B/A 2 , the relative strengths of the critical decay 
term (B) and Cole-Davidson term (A 2 ), by the follow- 
ing strategy. For fixed r, the location of the minimum 
in the function Lam" (w) occurs at a frequency io m m that 
decreases as B/A 2 increases. We then invoked the MCT 
result that io m in (unlike r a ) should be equivalent for all 
variables, and we adjusted B/A 2 to make to min fall close 
to the Lo m in found in the depolarized backscattering PC 
X"M spectra of Du et al [Du 1994]. For T = 190, 210, 
230, 250 and 270 K, the io m i n values of Du et al were 
found to agree quite well with the minima of the torn" (u) 
curves computed with Eq. (7) with B/A 2 = 0.175. With 



this fixed value of B/A 2 — 0.175, there are only two free 
parameters left in Eqs. (7) and (3): the CD. coupling 
constant A 2 , and the overall normalization constant Iq. 

With this severely constrained fit procedure, however, 
the resulting fits were not generally acceptable. One sur- 
prising aspect of the results was that for T < 170 K the 
linewidths of the fits became too broad; this problem can 
be resolved by decreasing B/A 2 . 

This indication that the strength of the critical contri- 
bution to the memory function (relative to the a peak) 
may weaken in the temperature range between Tc and 
Tq is also supported by the depolarized backscatter- 
ing spectra, shown in Fig. 7. At temperatures below 
T = 170 K, the minimum in loI{lo) is at lo < 1 GHz, so 
the a peak is completely out of the spectral window. The 
intensity in the Brillouin frequency region (~ 8 GHz), 
which is therefore entirely due to the "fast beta decay" , 
is seen to decrease rapidly with decreasing T, consistent 
with the apparent behavior of u>m(u>). This result is also 
presumably related to the MCT prediction of a cusp in 
the non-ergodicity parameter f q (T) at Tc below which 
the critical decay weakens with decreasing temperature. 
It also indicates that additional contributions to the fast 
relaxation from other processes (e.g. intramolecular vi- 
brations) are unlikely in PC since these would produce a 
larger value of B/A 2 . 
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FIG. 8. (HYBFIT4). Fits to PC Brillouin spectra 
Ijsoi 1 ^) with the hybrid model for m(u>) of Eq. (7). 
The free fitting parameters are lo, A 2 , and B/A 2 ; the 
values found from the fits are given in Table II. For 
T > 210 K, A 2 and B/A 2 were too strongly correlated 
to be determined by the fit, so B/A 2 was arbitrarily 
fixed at 4.0. 



E. Thermal diffusion contribution 



D. Complete hybrid model fits 

To provide a complete hybrid model analysis, we al- 
lowed the ratio B/A 2 to vary in the fit. At low tem- 
peratures, B and A 2 could be determined independently 
from the fit. At temperatures above 200 K, B and A 2 are 
strongly correlated. At 210 K an optimum fit was found 
with B/A 2 = 4.0, and we arbitrarily kept this value fixed 
for all temperatures above 210 K. The resulting set of fits 
is shown in Fig. 8; the fitting parameters B/A 2 and r 
are given in Table II together with the \ 2 values. These 
hybrid model fits have three free parameters, as in the 
CD. fits. However, they have t(T) constrained to follow 
Eq. (1) shown by the broken line in Fig. 1, automatically 
eliminating the t(T) disagreement problem of the CD. 
fits. They are nevertheless at least as good as the CD. 
fits. 

In their light scattering study of meta-toluidine, 
Aouadi et al [Pick 2000] also analyzed Iiso(w) spectra 
with the hybrid model, constraining t^ d to be propor- 
tional to the r a determined from depolarized backscatter- 
ing spectra. They found, as we did, that at low tempera- 
tures the strength of the critical decay decreases with 
decreasing T (see their Fig. 12). Although they also 
found that the critical decay strength decreases at high 
temperatures, resulting in an apparent maximum around 
250 K, this behavior may be a manifestation of the 
strong correlation between A 2 and B mentioned above. 



With a e ff replacing a in Eq. (7), we carried out ad- 
ditional fits to the Brillouin spectra for T — 150, 160, 
170, and 180 K. The fits are shown by the lower curves 
for each temperature in Fig. 9. (They are also the four 
low-T fits shown in Fig. 8 ) . Note that at low tempera- 
tures the experimental spectra have weak extra structure 
at low frequencies which is not present in the theoretical 
fits. 

The light scattering spectrum of density fluctuations 
Iisoi^), including thermal diffusion effects, can be ex- 
pressed as [Boon 1980], [Scopigno 2000] 



I(u>) = — Im [ui T — lo — iujjo — i 



LO 



»M]' 



(10) 



where lot = CtQ {Ct is the isothermal sound velocity) 
and m(ui) = m v (u>) + rriThi^) includes both the struc- 
tural relaxation contribution m v (io) [e.g. Eqs. (4) or (7)] 
and the thermal diffusion contribution 

m Th {t) = fr-l^e-*/™ 

where rrh — {Dt<I 2 )~ 1 ■ With the Laplace transform 
convention of Eq. (18) in the Appendix, 



m Th (u) 



z(7 - l)L0 2 TTh 
1 - lLOT Th 



(11) 



In the spectral window explored by 90° Brillouin scatter- 
ing (to > 0.5 GHz), u/TTh >> 1 since typically TTh ~ 



10 



15 ns. Then, in Eq. (10), —um T h{u) ~ ^5,(7 - 1). 
Adding this to the uj\ term in Eq. (10) then gives 
ll>\ + w|,(7 — 1) = cl>^7 which is just the square of the 
adiabatic sound frequency u)q. In that approximation 
Eq. (10) recovers the simpler equation [Eq. (3)] with 
m(u/) = m v (oj) only. 

While the quasielastic thermal diffusion mode is not 
usually included in the analysis of Brillouin scattering 
spectra, it is included in the analysis of inelastic x-ray 
scattering experiments since with larger q it becomes 
much broader [Scopigno 2000]. Also, at low tempera- 
tures, the high-frequency tail of the thermal diffusion 
mode can appear in the Brillouin spectrum if the window 
extends to sufficiently low frequencies. The four spectra 
shown in Fig. 9 have weak extra structure at low frequen- 
cies which may be due to the tail of the thermal diffusion 
mode. 

Since numerical data for 7 and Dj are not available 
for PC (to our knowledge), we estimated 7 = 1.5 and 
Trh = 16 ns, and reanalyzed the spectra of Fig. 9 using 
Eq. (10) with io\ = jlo 2 , m{u>) — m v (w) + m,Th{u) with 
rriThiw) given by Eq. (11), and m v (uj) by Eq. (7). The re- 
sulting fits, including thermal diffusion, are shown as the 
upper curves for each temperature in Fig. 9 for T = 150, 
160, and 170 K, showing that at these low temperatures, 
the thermal diffusion mode is apparently visible in the 
Brillouin spectrum at low frequencies. For temperatures 
above 170 K, the thermal diffusion contribution is not 
visible. 
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FIG. 9. (HYBFIT3). Low-temperature hybrid model 
fits using the a e ff values found in the fits shown in 
Fig. 7. For T = 150, 160, and 170 K, a second fit (the 
upper curve at each temperature) is included for each 
spectrum which includes the effect of thermal diffusion. 



F. Hybrid model with (3 — b 



Finally, we note that a slightly different strategy can be 
followed in carrying out the hybrid model analysis. With 
the exponent a in Eqs. (6) and (7) fixed (for T > Tq) at 
0.29, the region of the minimum in Lum"(uj) is predicted 
by MCT to obey Eq. (IV.5) with b = 0.50. However, 
there is a slight difference between this prediction and 
Eq. (7) because so far we have taken (3 ^ b. This differ- 
ence can be eliminated by arbitrarily fixing (3 = b = 0.50, 
as was done for the analysis of depolarized backscattering 
spectra of toulene by Wiedersich et al [Wiedersich 2000] . 

We therefore carried out another set of hybrid model 
fits with = b = 0.5. All other parameters were the same 
as in the fits of Fig. 8. In these fits, we kept a = 0.29 
and (3 = 0.50 for all temperatures and also kept both 
A 2 and B free. The resulting fits are comparable to 
those obtained with (3 = 0.68, demonstrating that the 
hybrid model fits are not very sensitive to the value of (3. 



At low temperatures, however, the low-frequency re- 
gions of these fits are not as good as the previous fits, 
illustrating the effect of the increase in a e ff at low T. 
The x 2 values obtained both with /3 = 0.68 and (3 = 0.5 
are listed in Table II together with the B/A 2 values, and 
are seen to be similar. Note that in the (3 = 0.5 fits where 
A 2 and B were both free fitting parameters, B/A 2 in- 
creases with increasing T, passes through a maximum at 
T ~ 230 K, and then decreases again at high T, similar to 
the results of Aouadi et al for meta-toluidine [Pick 2000] . 
As noted above, however, B and A 2 are strongly cor- 
related at high temperatures, and the high-T values of 
B/A 2 are therefore not significant. 

An advantage of choosing (3 = b is that there is one 
less free parameter in the fitting procedure. Also, the 
region of the susceptibility minimum in Eq. (9) is the 
same as in Eq. (7) if j3 = b. This makes the (3 = b hybrid 
model a good candidate for attempting to simultaneously 
describe all relaxing properties in a given material. 
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V. DATA ANALYSIS 3: EXTENDED 
SCHEMATIC MCT MODEL 

For the third and final analysis of the Brillouin spec- 
tra //5o(oj), we constructed the memory function m(u>) 
of Eq. (3) from the mode coupling theory rather than 
from parameterized empirical models, using the extended 
schematic MCT model of Gotze and Voigtmann men- 
tioned in the Introduction [Gotze 2000] . 



A. Mode coupling theory 

Mode coupling theory (MCT) begins with exact classi- 
cal equations of motion (the generalized Langevin equa- 
tion) for 4>q(t), the normalized autocorrelation functions 
of the density fluctuations p q (t) , derived from the Hamil- 
tonian with the Zwanzig-Mori projection-operator for- 
malism [Gotze 19921: 



4> q (t) + i/^,(t) + fi^,(t) + n l m S~ t')4>S')dt' = 

Jo 

(12) 

Equation (12) is equivalent to the generalized hydro- 
dynamics equation (3), as shown in the Appendix, but 
its parameters are determined by the Hamiltonian in the 
Zwanzig-Mori formalism. In particular, m q (t — t'), which 
is an empirical memory function in Eq. (3), is the auto- 
correlation function of the random forces in Eq. (12). A 
central accomplishment of MCT was the derivation of a 
closed mode-coupling approximation for m q (t) as 



i(t) = ^2 V ^ 9 1 ' q ~ qi ^<n (*)0<?-<?i (*) 



(13) 



<n 



in which the coupling coefficients V(q, q\, q — q\) are ex- 
pressed in terms of the static structure factors S q , S qi , 
and S, 



q-qi- 



If the intermolecular potential is known, 
Eqs. (12) and (13) can be solved self-consistcntly to 
find the density correlation functions and spectra [Fra- 
nosch 1997b]. 

In cases where the intermolecular potential is not 
known, simplified "schematic" MCT models can be used 
which include only a few correlators, with the remaining 
coupling constants considered as fitting parameters. The 
Sjogren model is a two-correlator schematic model. The 
"system" is represented by the system correlator <f>(t), 
while the variable being probed experimentally is repre- 
sented by the probe correlator 4> s (t). The two memory 
functions arc 



m{t)=v l (j){t)+V24> 2 {t) 



m s {t) = V s (j){t)(j) s {t) 



(14a) 



(14b) 



Equation (12), written separately for <j>(t) and 4> s (t) with 
the memory function of Eqs. (14), can then be solved 
sclf-consistently with v\, i>2 , and V s as fitting parameters 
(see [Singh 1998]). 

For temperatures near and below Tq of MCT, the 
"cage effect" dynamics embodied in Eq. (12) lead to the 
(usually) non-physical result of total structural arrest. 
In most glass-forming materials, this arrest is avoided 
by activated hopping processes which are included in the 
extended version of MCT (to a first approximation) via 
a temperature-dependent hopping parameter A. Equa- 
tion (12) is replaced by [Gotze 2000] 

4>(t) + (A + v)4>{t) + (fl 2 + A • v)<f>{t) 

+n 2 [ m(t - t')[<j>(t') + A • <f>(t')]dt' = (15) 

Jo 

(Note that the hopping parameter A in Eq. (15) is not 
related to the coupling constant A 2 of the Brillouin fits 
in, e.g., Eq. (7).) 

The extended Sjogren model of Gotze and Voigtmann 
consists of two Eqs. (15), one for <p and one for <p s , with 
the two memory functions of Eqs. (14). These arc to be 
solved self-consistently. 

The parameters of <j>(t) are: Q, v, A, v\, and i>2- They 
have already been globally optimized for PC by Gotze 
and Voigtmann for dielectric, neutron scattering and de- 
polarized light scattering data covering a much larger 
range of frequencies than that of our Brillouin scattering 
spectra. The parameters of 4> s {t) are: f2 s , u s , A s , and 
V s . (Note that varying the parameters of <j> s has no effect 
on <f>.) 

In this section, the Gotze- Voigtmann schematic MCT 
model will be used to compute the memory function 
m(oj) of Eq. (3), taken as m(u>) — A 2 tp s (uj), where A 2 
is the coupling constant. [We note that this is the same 
procedure employed by Ruffle et al [Ruffle 1999] who used 
the simpler Sjogren model)]. 



B. MCT fitting procedure 

Fitting the Brillouin spectra Iisoi^- 1 ) with the two- 
correlator MCT model described above is technically 
more difficult and therefore less straightforward than the 
two other fitting procedures. This is because m(u>) = 
4> s (u>) must first be obtained by solving and Fourier trans- 
forming the MCT equations [Eqs. (15)], and then interpo- 
lated to match the frequencies of the Iiso 0*0 spectrum. 

For the fits with the schematic MCT model, we began 
with the parameters obtained in the depolarized light 
scattering (DLS) fits of Gotze and Voigtmann, shown in 
their Figure 1 [Gotze 2000] . Since the correlators 4> s com- 
puted with the parameters they used gave very good fits 
to the depolarized light scattering spectra over nearly 
four decades in frequency, we tried to fit the Brillouin 
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scattering spectra using these same parameters. How- 
ever, the a peaks of the DLS spectra are at lower fre- 
quencies than the corresponding t^ a values, as seen in 
Fig. 1, where t® ls ~5x t^ a . We found that the po- 
sition of the a peak in uKp"(u}) is primarily determined 
by the coupling constant V s in Eq. (14b). We therefore 
treated V s as a free fitting parameter. All other MCT 
parameters were kept at their DLS values. We solved 
the MCT equations [Eqs. (15)] for selected values of V s 
and the resulting cf> s (u>), after multiplication by the (ad- 
justable) coupling constant A 2 , was taken as the memory 
function m(uj) in Eq. (3). A nonlinear least-squares anal- 
ysis was then carried out in which Iq and A 2 were the free 
parameters. This process was repeated until the fit was 
optimized. (Note that in fitting their Brillouin spectra 
of NaLiPC>3 with the idealized two-correlator schematic 
Sjogren model, Ruffle et al also treated V s as a free pa- 
rameter while the system parameters were fixed from 
neutron scattering results [Ruffle et al 1999].) 

The MCT fits for T = 160, 180, 200, 220, 240, 260, 
280, 300, 320, and 350 K are shown in Fig. 10. (For 
T > 300 K, the MCT parameters other than V s were ob- 
tained by extrapolation from lower temperatures.) The 
values V s and x 2 for these fits are given in Table II (note 
that the V s values for the Brillouin fits are systemati- 
cally smaller than the depolarized light scattering val- 
ues V S DLS ). The MCT fits are comparable to the hybrid 
fits, and have lower x 2 values at some temperatures. For 
the lowest temperatures, the low-frequency region of the 
MCT fits is poor, reflecting the known limitations of the 
model for low frequencies at T < T c [Gotze 2000b]. At 
temperatures above 280 K, the fits were not sensitive to 
the value of V s for V s < 3.5. This is because the a peak in 
u>m"(uj) is then above the Brillouin line, and the shape of 
the memory function in the Brillouin range (1 to 15 GHz) 
is nearly independent of V s . 

For the temperatures 180, 200, 220, 240, and 260 K, the 
position of the a peak in wm"(w) was sufficiently well- 
determined by the fits to estimate the MCT a-relaxation 
time r a — 1/wq from the frequency of the a peak. The 
values found are listed in Table II and are also shown in 
Fig. 1 by the * symbols. Note that the MCT values follow 
the same qualitative temperature dependence as the di- 
electric data, and are similar to (though somewhat lower 
than) the values we assumed for the hybrid fits. Thus, 
the MCT approach, in contrast to the other approaches 
followed in the preceding two sections, is able to provide 
reasonable estimates of t^ a purely on the basis of fits to 
the Brillouin spectra. 

These results show that this extended schematic MCT 
model is able to simultaneously describe the spectra and 
also give reasonable fits for r a (T), thus answering the 
question that originally motivated this study affirma- 
tively. 



|MCTIil3.pll1 [9 Jul) 01 1| 







350 
















1 o 3 - 










160 










1 D 4 -. 


lai# 






\^isv 












1 o 3 -. 


^|pp«*«=— ■ 


2-^ 


^ 






^Sk. 












■^ ■- #*~v 


yga 


■«$>& 










J 5 .: 




10' - 1 




















} 2 4 




6 


8 


10 


12 


14 




16 



FIG. 10. (MCTFIT3). Fits of PC I ISO (w) spectra with 
the extended schematic MCT model of Eqs. (15) and 
(14) for T = 350, 320, 300, 280, 260, 240, 220, 200, 180, 
and 160 K. In the fits, j (T) = T/1000 and w (T) is 
fixed following Eq. (5). All MCT parameters are those 
found by Gotze and Voigtmann from their fits to DLS 
spectra except for V s which has been adjusted to op- 
timize the fits. For T > 300 K, the parameters other 
than V s were determined by extrapolation from lower 
temperatures. 



VI. DISCUSSION AND CONCLUSIONS 

As noted over 30 years ago by Montrose et al, the de- 
tailed shape of the Iiso^) Brilliouin spectrum is rel- 
atively insensitive to the exact nature of the memory 
function [Montrose 1968]. This is due, as noted ear- 
lier, to the limited frequency range of Brillouin scatter- 
ing data. If wo is fixed (from ultrasonic data), the ratio 
m" '(u> b) I 'tu' '(to b) must have the "correct" value to obtain 
reasonable fits, but there is still enormous flexibility in 
constructing m(u>). If one wishes to extract meaning- 
ful information from the fits, then additional constraints 
must be imposed. 

First, the a relaxation time should be reasonable. Fits 
using a-relaxation-only models, e.g. the Cole-Davidson 
function of Eq. (4), produce acceptable fits but give un- 
reasonably short relaxation times at low temperatures. 
At T ~ T G , one finds typically t^ d (T g ) ~ 10" 7 sec 
while other experimental techniques generally yield 
T a{Tc) ~ 10 2 sec. This result demonstrates the fact, 
already shown many times, that a-only models of struc- 
tural relaxation are incomplete. 

A great deal of research during the past 15 years has 
shown that structural relaxation occurs in two steps, 
starting with a fast decay towards a plateau that precedes 
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the final a-relaxation process. Since the memory func- 
tion m(oj) for longitudinal viscosity should also have this 
two-step form, the memory function used in the analysis 
of Brillouin scattering spectra should include both the a- 
rclaxation and the fast-relaxation parts of m(u>). If r Q is 
constrained to follow the temperature dependence found 
in other measurements, such as dielectric or depolarized 
light scattering spectroscopy, then the memory function 
must be extended to higher frequencies to provide suf- 
ficient damping once the a peak in Lum"(Lu) has moved 
below the Brillouin lines. This extension to include fast 
relaxation processes can be implemented many ways, but 
most approaches require the introduction of a fast j3 re- 
laxation time T/3 that does not correspond to the results 
of other experiments. 

In practice, there may also be other sources of fast re- 
laxation - such as intramolecular vibrations - that influ- 
ence the Brillouin spectra, but the high-frequency part of 
the structural relaxation dynamics must still be present. 

Because the mode coupling theory has been successful 
in analyzing a wide range of experimental data, we have 
used an empirical model for m(w) which is qualitatively 
consistent with MCT. The hybrid model, which mimics 
the two-step relaxation scenario of MCT by combining 
the CD. function for a relaxation with the critical decay 
law (lu 11 ) has been shown to provide generally excellent 
fits to the Brillouin scattering spectra. 

Furthermore, we found that fixing [3 = b = 0.5 pro- 
duces fits of equally good quality as with (3 = 0.68. This 
reduction in the number of free fitting parameters rep- 
resents an attractive variant of the hybrid model which 
has been shown, for the case of toluene, to provide good 
fits to the depolarized backscattering spectra [Wiedersich 
2000] . It will be interesting to explore the ability of this 
reduced hybrid model to simultaneously fit both depo- 
larized backscattering spectra and Brillouin spectra for 
any given material. 

We have also explored the applicability of the two- 
correlator schematic MCT model to the analysis of Bril- 
louin spectra, a procedure first employed by Ruffle et al 
[Ruffle 1999] who used the Sjogren model without hop- 
ping. We used the new extended model of Gotze and 
Voigtmann [Gotze 2000] with which they fit depolarized 
light scattering, neutron scattering, and dielectric data 
for propylene carbonate with common parameters for 
the system correlator <p(t) for all the experiments. For 
the probe correlator <f> s (t), we kept the parameters Q s , 
v s , and A s they obtained from fits to depolarized light- 
scattering spectra fixed, but we varied V s to optimize 
the fits to the Brillouin spectra, taking m{uj) = A 2 (f> s (<jj) 
where A 2 is a free parameter. The resulting V s values 
were systematically smaller [and the resulting frequency 
of the a peak in com" (to) higher] than those found from 
fits to the depolarized light scattering spectra, consistent 
with the result that t^ a is smaller than r a of depolarized 
light-scattering spectra as shown in Fig. 1. We then used 



the positions of the a peaks in the MCT memory func- 
tion to estimate the relaxation time T^f CT = 1/ui^f . 
As seen in Fig. 1, the resulting t^ ct values follow the 
temperature dependence expected on the basis of the di- 
electric data of Schneider et al. We emphasize that this 
MCT analysis is the only one of the three employed in 
this study that can provide a meaningful estimate of r Q . 

The recent Brillouin scattering study of OTP by 
Monaco et al [Monaco 1999, Monaco 2000, Monaco 2001, 
Mossa 2001] has suggested a different origin for the 
fast relaxation process observed in Brillouin scattering 
spectra: coupling of longitudinal acoustic waves to in- 
tramolecular degrees of freedom as originally suggested 
by Mountain [Mountain 1968]. Even in OTP, however, 
where this mechanism appears to be important, some 
part of the fast relaxation must be associated with struc- 
tural relaxation. It is therefore an open question, that 
should be carefully addressed, how much of the fast re- 
laxation originates in each of these two mechanisms, in 
OTP as well as in other materials - including PC. 

Finally, we return to the question of the ability 
of Brillouin scattering spectra to extract meaning- 
ful information about the structural relaxation pro- 
cess. In Fig. 11 we show ujm"{uj) (solid lines) 
and ujm'(uj) (broken lines) obtained from the MCT 
fits (top) and the hybrid fits with j3 = b (bottom). 
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FIG. 11. (MEMPLOT200). wm"(w) (solid lines) and 
com'(uj) (broken lines) from the fits at 200 K. Top: 
MCT fit; bottom: hybrid fit with (i = b = 0.5 and 
a = 0.29. The fit range (1-15 GHz) is indicated in both 
plots by a thick horizontal line. 
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In the range of the Brillouin spectra (1-15 GHz indi- 
cated by the horizontal solid lines), the shapes of these 
curves are remarkably close, even though they were ob- 
tained with very different procedures. This result sug- 
gests that a reasonable phcnomcnological model like the 
hybrid model, with the constraint that r a should have a 
temperature dependence consistent with that found by 
other techniques, can provide a credible representation 
of the memory function m(u>) characterizing the longitu- 
dinal viscosity. 
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VII. APPENDIX: SOME MATHEMATICAL 
DETAILS 

The normalized density fluctuation correlation func- 
tion 



<j>(t) =< Pq (0) Pq (t) > / < \ Pq (0)\ 2 > 



(16) 



obeys the generalized Langevin (or generalized oscillator) 
equation 



4>(t) + 7o^(i) + ul4>{t) + m(t- t')<j)(t')dt' = 



(17) 



We use the Laplace transform convention usually fol- 
lowed in the mode coupling theory literature 



F{z) = i / e lzt F{t)dt 
Jo 

with which we obtain 

-[z + m'(z)} -i[yo + m"{z)] 



(z) 



[z 2 — luq + zm'(z)] + i[z"fo + zm"(z)] 



(18) 



(19) 



where m(z) — ra'(z) + im"(z). 

The power spectrum of the fluctuations in <fi(t) is given 
by 

S(w) = Im[4>(z)] z = u (20) 

and the light scattering spectrum I(u>) oc S(u>) is then 

7 [7o+m"(w)] 



/(«) = 



[lo 2 — u>q + uim'(uj)] 2 + [lo^q + wm"(w)] 2 



(21) 



which is Eq. (3) in the text. (Note that without the i in 
Eq. (18), to' and to" would be exchanged, a convention 
that is frequently employed in the Brillouin scattering 
literature.) 

In the critical region preceding the a relaxation m(t) 
can be represented by 



m cri t (t) =t e 



— +- a p -t/ T c 



(22) 



wto(w) = A 2 [(l - iurY !i - 1] + iwB x 

= A 2 [(l - iojT)-P - 1] + iiuB x T(l - a)T l - a (l - iLur) 



when the e~ t ' Ta cutoff factor prevents m C rit{t) from 
adding to m a (t) beyond the center of the a relaxation 
region. With Eq. (18) and m(u>) = [m(z)] z=w we have 

m crit (uj)=i e luJt r a e' t/T dt (23) 

Jo 

from which [Gradshteyn 1965, p. 317] 

m crlt {u) = iT(l-a)(--iLj) a - 1 (24) 

T 

For the hybrid model, m{ui) = mcD{u) + m cr i t (uj), so 
T(l - a)[r- 1 - iu] a - 1 



a-l 



(25) 



which is Eq. (7) in the text. (For a = 0.29, T(l - a) 
1.282.) 
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TABLE I. Properties of Propylene Carbonate 



Property 


Value 


Reference 


Formula 


C4H6O3 




Molecular weight 


f02.09 




Tboil 


513 K 




Tm 


218 K 




T g 


160 K 




T c (MCT) 


185 ± 5 K 


Gotze a 


Density p(T) 


1.541 - 1.148 x 1(T 3 T (K)(g/cm 3 ) 


Sim b 


Refractive index Hd(T) 


1.5314-3.752 x 1(T 4 T (K) 


Sim 


Sound velocity C (T) 


(2.5075 x 10 5 - 361.2 x T) cm/sec 


Du c 


Pk (Kohlrausch stretching coefficient) 


0.77 ±0.05 


Du 


Pcd (Cole-Davidson stretching coefficient) 


0.68 


Du 


Tc (MCT crossover temperature) 


187±5K 


Du, Sch d 




182 


Wut c 


A (MCT exponent parameter) 


0.78 ±0.5 


Du 




0.72 


Wut 


m (fragility index) 


104 


Bohm f 



a (Gotze) W. G6tze and Th. Volgtmann, Phys. Rev. E 61, 4133 (2000). 

b (Sim) L. Simeral and R.L. Amey, J. Phys. Chem. 74, 1443 (1970). 

c (Du) W.M. Du, G. Li, H.Z. Cummins, M. Fuchs, J. Toulouse, and LA. Knauss, Phys. Rev. E 49, 2192 (1994). 

d (Sch) U. Schneider, P. Lunkenheimer, R. Brand, and A. Loidl, Phys. Rev. E 59, 6924 (1999). 

c (Wut) J. Wuttkc, M. Ohl, M. Goldhammer, S. Roth, U. Schneider, R. Lunkenheimer, R. Kahn, B. Ruffle, R. Lechner, and 

M.A. Berg (submitted to Phys. Rev. E). 

f (Bohm) R. Bohmer, K.L. Ngai, C.A. Angell, and D.J. Plazek, J. Chem. Phys. 99, 4201 (1993). 



TABLE II. 



T(K) 


CD. (free r) 




Hybrid (fixed r) 


Schematic MCT | 






13 = 0.68 


/3 = fe = 0.5 


t (ns) 


* s 


yBH 


x 2 




t (ns) 


A" (GHz^) 


x 2 


r (ns) 


B/A' 


x 2 


B/A* 


x 2 










140 


790 


36.5 


27.2 


1.26 x 10^ 


















150 


1220 


40.2 


27.8 


1.42 x 10 1 * 


0.01 


27.4 


0.002 


28.6 




94 






160 


31.8 


37.0 


22.1 


2.49 x 10 B 


0.04 


19.5 


0.034 


22.7 




80 


30 


25.2 


165 


14.2 


34.3 


23.3 


1.96 x 10° 


















170 


6.11 


31.8 


41.3 


4.28 x 10" 


0.11 


12.5 


0.13 


15.9 




60 






180 


1.34 


27.5 


66.3 


2.11 x 10^ 


0.53 


8.4 


0.32 


20.3 


2.9 x 10 1 


50 


13 


31.4 


185 


0.763 


25.3 


84.6 


3.40 x 10 1 


















195 


0.236 


22.0 


60.4 


2.51 


1.50 


47.3 














200 


0.148 


20.6 


44.2 


9.75 x 10 1 




35.2 


1.31 


23.3 


2.4 x 10 _I 


30 


7.6 


5.85 


205 


1.03 x 10 _1 


19.5 


28.8 


4.47 x 10 _i 


















210 


7.28 x 10"^ 


18.4 


16.4 


2.33 x 10 _i 


4.00 


10.1 


2.72 


8.3 




25 






215 


5.41 x 10^ 


17.6 


17.8 


1.34 x 10 _i 


















220 


4.15 x 10"^ 


17.0 


8.99 


8.35 x 10^ 




19.2 


18.1 


5.75 


3.0 x 10^ 


23 


5.4 


5.73 


225 


3.25 x 10"^ 


16.9 


9.43 


5.55 x 10^ 


4.00 


19.5 














230 


2.53 x 10"^ 


16.7 


5.21 


3.89 x 10^ 




20.3 


36.0 


12.1 




21 






235 


2.13 x 10"^ 


17.0 


5.35 


2.85 x 10~" 


















240 


1.75 x 10^ 


17.0 


11.9 


2.16 x 10~^ 


4.00 


14.2 


24.8 


13.2 


9.9 x 10~' 1 


18 


4.7 


11.9 


250 


1.28 x 10^ 


17.0 


25.7 


1.36 x 10"" 
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260 


9.96 x W-* 


17.0 


50.0 


9.29 x lO - ' 1 


4.00 


41.6 


21.2 


41.3 


4.7 x 10~' 1 


15 


4.3 


39.9 


270 


8.17 x 10~ a 


17.0 


47.4 


6.80 x 10- 3 












13 






280 


6.61 x 10~ rt 


17.0 


66.3 


5.24 x 10~ J 


4.00 


54.3 


22.0 


54.0 




12 


<3.5 


51.9 


290 


5.58 x W-* 


17.0 


60.5 


4.19 x 10~ J 


















300 


4.57 x 10~ rt 


17.0 


60.7 


3.46 x 10~ J 


4.00 


53.5 


0.87 


53.5 






< 2.5 




310 


3.90 x 10~ a 


17.0 


43.4 


2.93 x 10~ J 


















320 


3.38 x W-* 


17.0 


49.2 


2.53 x 10~ J 


4.00 


44.3 


2.34 


44.3 






< 2.5 


40.5 


330 


2.76 x 10~ a 


17.0 


38.1 


2.22 x 10^ 


















340 


2.26 x 10~ rt 


17.0 


26.3 


1.98 x 10~ J 


















350 


1.80 x 10~ a 


17.0 


18.2 


1.79 x 10~ J 


4.00 


18.1 


1.90 


18.1 






< 2.5 


18 



17 



